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We calculate the zero temperature (T = 0) phase diagram of a polarized two-component Fermi gas 
in an array of weakly-coupled parallel one-dimensional (ID) "tubes" produced by a two-dimensional 
optical lattice. Increasing the lattice strength drives a crossover from three-dimensional (3D) to ID 
behavior, stabilizing the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) modulated superfluid phase. We 
argue that the most promising regime for observing the FFLO phase is in the quasi- ID regime, where 
the atomic motion is largely ID but there is weak tunneling in the other directions that stabilizes long 
range order. In the FFLO phase, we describe a phase transition where the quasiparticle spectrum 
changes from gapless near the 3D regime to gapped in quasi-lD. 



Recent experiments on ultracold ^Li have probed po- 
larized, two-component Fermi gases as a function of in- 
teratomic interaction strength and spin population im- 
balance [H . These experiments have focussed on the uni- 
tarity regime, where the s-wave scattering length is large 
and the scattering properties are universal. At low tem- 
peratures, they have seen phase separation between an 
unpolarized or weakly polarized superfluid phase and a 
highly polarized normal phase. Future experiments hope 
to observe the elusive FFLO spatially-modulated super- 
fluid phase, first predicted to occur in magnetized super- 
conductors over 40 years ago [3, Q i and having realiza- 
tions in other systems ranging from heavy-fermion su- 
perconductors Q to quark matter While the FFLO 
phase is expected to exist in trapped 3D gases for small 
polarizations and weak attractive interactions it is 
predicted to only occupy a small part of the T = phase 
diagram Q, and this region is only further diminished 
with increasing temperature [8|]. Here we theoretically 
explore how a two-dimensional (2D) optical lattice can 
enlarge the region of stability of the FFLO phase, paving 
the way for its observation. 

Although a 3D simple cubic optical lattice may also 
enhance the FFLO phase 0, we argue that a 2D opti- 
cal lattice, which permits free motion in one direction, 
is more effective. Analogous to charge density wave in- 
stabilities, the instability of the normal state to FFLO is 
due to a Fermi surface "nesting" , which is enhanced in 
ID. By increasing the intensity of the optical lattice, one 
can continuously tune the single atom dispersion from 
3D to ID, a scenario which is readily achieved experi- 
mentally [10|. 

As revealed by the Bethe ansatz 11, [13], the exact 
T = phase diagram of the ID polarized Fermi gas dis- 
plays four phases: unpolarized superfluid (SF), FFLO, 
fully-polarized normal, and vacuum, characterized by the 
densities of the two species and by algebraic order. Unlike 
3D, in ID all of the phase transitions are continuous and 
the FFLO phase occurs at all non-zero partial polariza- 
tions for any strength of the attractive interaction [l^ . 



Furthermore, at the phase boundary between polarized 
and unpolarized phases, the SF phase has the lower den- 
sity in ID, which is opposite to the situation in 3D. In a 
trap, the spatial structure is therefore inverted: in ID one 
has a central FFLO region surrounded by SF, for small 
polarizations. During the crossover from ID to 3D, we 
find a regime where the phase sequence moving from the 
center to the edge of the trap is the quasi-lD FFLO, then 
a shell of SF, then the more 3D-like polarized normal. 

We consider a gas of fermionic atoms in two dif- 
ferent hyperfine states (labeled by a =T,i= ±1) 
confined by a smooth trapping potential V{y) and 
a square optical lattice potential such as Vl{y) = 
—U{y) [cos(27ra;/6) -I- cos(27r?//&)] , with lattice spacing b 
and local depth ?7(r), which breaks the cloud into an 
array of tubes aligned along the z-direction. In the lo- 
cal density approximation (LDA), the properties at loca- 
tion r depend on U (r) and the local chemical potentials 
/io-(r) = /Ict — V{v) in the same way that a spatially uni- 
form system does. Within LDA, which should be valid 
for a wide enough trap, the spatially varying pattern of 
phases in a trap can be read ofl from the phase diagram 
of the uniform system by tracing the spatial variation 
of /i = -I- Hi)/ 2 and U , while holding the difference 
h = (/i| — iJ.i)/2 flxed. In the special case of uniform U, 
only fi varies in space. 

To produce the uniform phase diagram we study the 
untrapped system in a uniform lattice with x Ny 
tubes, each of length in the z-dircction. For a suf- 
ficiently low density and strong enough lattice, the xy 
motion is well approximated by a one-band tight-binding 
model with single-atom dispersion: 
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2t[2 — cos{kxb) — cos(A:j,6)] 
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where t is the hopping (related to U and b as in 
m is the atomic mass, and we use h — 1. The 
Brillouin zone of the xy motion is |fca;|,|fcy| < 7r/6, 
while kz is unconstrained. For energies well above 
the xy bandwidth ek ^ 8t, the dispersion is ID-like. 
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For low energies ek <C t, the dispersion is 3D-like, 
and can be made isotropic if we rescale the momenta 
{k^, ky, kz} {b^/2ik.j;,b^/2iky,kz/y/rn}. This single- 
band, tight-binding regime is accessed experimentally by 
working in a regime with t^epa <SC \JU jmlP- ^ where 
£_Fcr = {unch'^Y /2m is the ID Fermi energy for each 
species of density (corresponding to the ID density 
ricrb^ per tube in the optical lattice). 

Since the ^Li experiments use highly dilute gases with 
a wide Feshbach resonance, the interactions can be mod- 
eled by a contact interaction, giving a Hamiltonian, 

H - A*T^T - Mi^i " H 51 ("^k - M^) cLck^ 

k <T=Ta 

^ X! CkT^'iCk'+qiCk-qT • (2) 



LzN^Ny 



k,k',q 



Solving the 3D scattering problem in a single harmonic 
tube of transverse size aj_ = -^/I/tocj^ , one finds 14 1: 
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where a is the 3D s-wave scattering length, arising from 
the short-range interatomic potential. Thus, we have an 
attractive ID interaction when aj_/a < C ~ 1.4603/V2. 
Defining the ID two-body binding energy es — g'^'m/A, 
we can fully parameterize the T = phase diagram of 
the uniform system with three dimensionless quantities: 
t/sB, IJ^/sb and h/eB |15l |. 

We calculate the T = phase diagram within mean- 
field theory, which captures most of the qualitative fea- 
tures of the phase diagram as we move between the ID 
and 3D regimes. However, we know by comparison to the 
exact solution that this mean-field approximation does 
miss some features of the ID limit, as we note below. 

We begin at /i = where there are only two phases: 
SF and the vacuum (see Fig. [T]). In the ID limit (small 
t) there is a two-atom bound state with binding energy 
Eb, where clearly Eb — £b when t = 0. These bosonic 
pairs enter the system and form a Bose-Einstein con- 
densate (BEC) for /i > —Eb/2. Increasing the density 
further brings the system through a density-driven BEC- 
BCS crossover, similar to excitonic systems .16j, where 
h/Eb ^ 1 defines the weak-coupling BCS limit. Mak- 
ing the system more 3D, by increasing reduces both 
Eb and the BEC regime. For I/eb > 0.2066 there is no 
two-atom bound state and thus only the BCS regime. 

At finite h we must also consider the FFLO superfluid 
phase, where the Cooper pairs condense with nonuniform 
pairing order parameter (gap) 
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FIG. 1: (Color online) Phase diagram for h — 0. For t/es 
below the filled circle, there is a two-atom bound state, and 
the resulting bosonic pairs enter the system as a BEC as /i is 
increased through the solid line (given hy ^ — —Eb/2). For 
t/sB above the filled circle we are always in the BCS regime. 



where f{x) is the zero-temperature Fermi function, k_L 
is the momentum in the xy plane, / labels the quasipar- 
ticle modes, and the energies/coherence factors obey the 
Bogohubov-de Gennes (BdG) equations [l3| 
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with horr = 



■'-Ot = + 2t[2 - cos{k^b) - cos{kyb)] + 

gn^„{z) — II. The densities in the Hartree term 
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n^iz) = (l/N^N, 
grand potential is [1 
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where the sum includes both positive and negative en- 
ergy eigenvalues. The simplest ansatz for the FFLO 
phase is Fulde and Ferrell's one-plane wave state A(z) = 
Here, the energy eigenvalues £'(,ki reduce 

to Ski - ±7(^i+d)V4 + A2 + (^1 - ei)/2, with 
— i]i±qz/2— IJ'+gn-a- and we can then minimize Eq. ^ 
directly. This state is a good approximation in the limit 
App 0. Indeed, one can show that the second-order 
transition to the normal phase occurs at single wave vec- 
tor q when [l^ 
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The locus of this transition is illustrated in Fig. [2l 

Larkin and Ovchinnikov 0] showed that the energy is 
lower if the Cooper pairs condense in a standing wave, 
with A(z) = Alo cos{qz) when the gap is small. More 
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FIG. 2: (Color online) Slice of the mean-field phase diagram 
taken at t/sB = 0.08. The phases shown include the unpo- 
larized superfluid (SF), partially-polarized normal (N), and 
fully-polarized normal (NP). The FFLO phase is divided into 
gapped 'commensurate' (C) and ungapped 'incommensurate' 
(IC) phases. The filled circle marks the tricritical point; near 
it, but not visible here is a tiny region of SFm magnetized 
superfluid phase, a remnant of the 3D BEC regime. The SF- 
NP and SF-N transitions are first-order for h/eb above the 
tricritical point, along the solid heavy line. The SF-FFLO 
transition (solid line) is estimated from the domain wall calcu- 
lation. The transition from FFLO to normal (dotted-dashed 
line) is assumed to be second-order. The large circle marks 
the region of FFLO where A/ef is largest, so the phase is 
likely to be most robust to T > here. The dashed line near 
the SF-FFLO transition shows where the wave vector of the 
FFLO state is stationary as a function of /j.: dq/dn — (this 
is calculated using the FF approximation). 



generally, A(z) is a real periodic function of z. When the 
coherence length ^ is small compared to l/q, this state 
consists of well-separated domain walls between domains 
where A is alternately positive and negative. The polar- 
ized cores of these domain walls result from occupying 
the spin- up Andreev bound states on each wall 21 1 . 

We calculate the energy of a single domain wall by it- 
erating to self-consistency Eq. ([S]) in a finite box with pe- 
riodic boundary conditions, beginning with a trial A.{z) 
containing two domain walls whose separation is large 
compared to the coherence length. If the domain walls 
interact repulsively, the SF to FFLO transition is con- 
tinuous and lies where this domain wall energy vanishes; 
otherwise this condition marks the spinodal of a first- 
order transition (likely to be near the true phase bound- 
ary). Within mean field theory the transition is contin- 



22l |. and has been argued to be so in 3D 



uous in ID 

pll 0, in weak coupling the critical fields are re- 
spectively h = (2/7r)Ao(= 0.64Ao), and 0.67Ao, where 
Aq is the gap in the SF phase. We are unaware of a 
strong coupling 3D calculation of the sign of the domain 
wall interaction. 



Fig. [2] shows a representative slice of the mean-field 
phase diagram at fixed t/eB = 0.08 (if one can neglect 
the spatial variation of and t, this slice corresponds 
to a fixed optical lattice intensity). Near the vacuum at 
small filling (low /x) is the 3D BEC regime, including a 
very small region of the SFm magnetized superfiuid phase 
where the excess fcrmions form a Fermi liquid within the 
BEC. As /i and thus the filling is increased, the system 
crosses over towards ID. Here, the FFLO phase appears 
and occupies a large portion of the phase diagram [25| . 
Both the SF and FFLO phases become re-entrant: in the 
ID regime the FFLO phase is at a higher /i and thus a 
higher density than SF, while in the 3D regime this den- 
sity relation is reversed. Thus, we see that the "inverted" 
phase separation in ID trapped gases is connected to 
the standard phase separation of 3D via an intermedi- 
ate pattern of phases where SF forms a shell surrounded 
by polarized phases. As t/es is further reduced, the 3D 
regime becomes smaller, with the re-entrance of the SF 
phase moving to lower ^, while the FFLO phase grows 
and the sliver of N phase between FFLO and NP is di- 
minished. In the limit t = this phase diagram matches 
fairly well to that obtained from the exact solution in 
ID (e.g.. Fig. 1 of Ref. The main feature that the 

mean-field approximation misses at t = is the multi- 
critical point where the four phases, SF, FFLO, NP and 
vacuum, all meet at /i = — /i = £b/2. In mean-field the- 
ory, the FFLO phase never extends all the way down 
to zero density; instead it is preempted by a first-order 
SF-to-NP transition. 

A new T = phase transition occurs within the FFLO 
phase as one moves from 3D to ID by increasing the 
intensity of the 2D optical lattice. In 3D the FFLO state 
has a Fermi surface, and is therefore gapless. In ID the 
spectrum of BdG quasiparticles is fully gapped in the 
FFLO state. The gapped, commensurate FFLO state 
(FFLO-C) contains exactly one excess spin-up atom per 
ID tube per domain wall. This commensurability means 
that q = 7r6^(ri| — «x), while, by contrast, the number of 
excess up spins in the ungapped, incommensurate FFLO 
state (FFLO-IC) is not constrained. 

The transition between FFLO-C and FFLO-IC can 
be understood from the band structure of the Andreev 
bound states on the domain walls. In FFLO-C the chem- 
ical potential lies in a gap in the quasiparticle spectrum. 
Thus, the superfiuid FFLO-C phase is a band insulator 
for the relative motion of the unpaired atoms and the 
condensate of pairs. As the optical lattice intensity is 
decreased, the 3D bands broaden and may overlap the 
chemical potential, opening up a Fermi surface. We ap- 
proximate the IC-C transition within the FF ansatz by 
examining the fc^ > half of the Fermi surface to see if 
it is fully gapped. In the limit ^/t ^ 1, the transition 
occurs when A ~ 8th/ ij. 

We now address the question of what are the best 
conditions for experimentally producing, detecting and 
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studying the FFLO phase. Ideally one might use in situ 
imaging to directly observe the spatial density and mag- 
netization modulations in this phase. In a trapped 3D 
gas, the modulated superfluid will occupy a hard to de- 
tect thin shell. The thinness of this shell results from the 
small range of fi over which the FFLO phase is stable 0] . 
Even approaches which produce an enlarged FFLO re- 
gion in density space [21| share this feature. Moreover, 
imaging the modulations will be complicated by their 3D 
nature (e.g., they may form an onion- like pattern). The 
ID limit also has problems. Although in ID the FFLO 
phase occurs in a large region of the T = phase dia- 
gram, it has no true long-range order, only power-law cor- 
relations. Furthermore, the transition temperature (Tc) 
of this ID superfluid is zero. Also, for t = (strictly 
ID) we have an array of independent parallel ID clouds 
whose density modulations will be out of phase with one 
another, reducing the observed signal. 

Given these concerns, we believe that the optimal con- 
ditions for observing FFLO are likely to be in the quasi- 
ID regime, where the 2D optical lattice is at an interme- 
diate intensity that is strong enough to make the Fermi 
surface ID-like (and hence enhance the instability to- 
wards FFLO), but weak enough that the atoms are still 
able to hop between the parallel tubes and thus intro- 
duce strong inter-tube correlations in the optical lattice. 
The resulting 3D long-range order can then survive to 
nonzero temperature [2J|. 

Although we have only performed aT ~ calculation, 
we can crudely estimate Tc from the size of the gap A. 
For small A, superfluid phases have Tc cx A, but when 
A approaches ep, Tc saturates. Thus the observability 
of a superfluid phase such as FFLO is enhanced if the 
gap is increased to of order sp, but there is not likely 
to be an advantage to increasing the gap to much larger 
values. In 3D the maximum Tc/sp of the SF occurs on 
the BEG side of the Feshbach resonance, well away from 
the FFLO phase [i^. However, as we move towards ID 
in our phase diagram, the FFLO phase extends more and 
more into the regime of strong pairing where the gap is 
of order ep, and thus we expect a large Tc- For a given 
t/sB, we find that the gap in the FFLO phase is the 
largest fraction of ep at the SF-FFLO phase boundary 
near its point of re-entrance, where h on the SF-FFLO 
boundary reaches its maximum (see Fig. 2). We also 
find that within mean-field theory this fraction A/ep in- 
creases as we reduce the hopping t. At sufficiently low 
t the system crosses over from quasi-lD to ID and our 
mean-field theory becomes unreliable. In the ID limit, 
Tc must vanish, so the maximum value of Tc/ep within 
the FFLO phase must occur in the quasi-lD regime at 
some small but nonzero hopping t. 

Another consideration that may complicate the detec- 
tion of the FFLO phase within a trap is the fact that 
fi varies spatially both within and between tubes. This 
means that the local wavenumber q of the modulation 



will vary through the cloud, making the modulations 
more difficult to detect. However, this variation can be 
minimized if one works near a point where dq/dfj, = 0. 
We find such points do exist in the quasi-lD regime (see 
Fig. (21); in 3D, dq/dfi is always negative so such points do 
not exist. Note that in the ID limit {t = 0) there is even 
a point in the exact phase diagram near strong coupling 
where d'^q/dfi = dq/d^ = that should be a real "sweet 
spot" for having a uniform q over a fairly large fraction 
of a trap, and that this feature should survive to small t. 
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